Construction of Equientangled Bases in Arbitrary Dimensions via 
Quadratic Gauss Sums and Graph States 
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Recently Karimipour and Memarzadeh [Phys. Rev. A 73, 012329 (2006)] studied the problem of 
finding a family of orthonormal bases in a bipartite space each of dimension D with the following 
properties: (i) The family continuously interpolates between a product basis and a maximally 
entangled basis as some parameter t is varied, and (ii) for a fixed t, all basis states have the same 
amount of entanglement. The authors derived a necessary condition and provided explicit solutions 
for D sC 5 but the existence of a solution for arbitrary dimensions remained an open problem. 
We prove that such families exist in arbitrary dimensions by providing two simple solutions, one 
employing the properties of quadratic Gauss sums and the other using graph states. The latter can 
be generalized to more than two parties - we can construct equally entangled families that vary 
continuously between a product basis and a graph state basis. 

PACS numbers: 03.67.Mn, 03.67.Hk 
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We present two different solutions to the problem 
posed by Karimipour and Memarzadeh in [1] of con- 
structing an orthonormal basis of two qudits with the 
following properties: (i) The basis continuously changes 
from a product basis (every basis state is a product state) 
to a maximally entangled basis (every basis state is max- 
imally entangled), by varying some parameter t, and (ii) 
for a fixed t, all basis states are equally entangled. As 
mentioned in [l[ , such a family of bases may find applica- 
tions in various quantum information protocols including 
quantum cryptography, optimal Bell tests, investigation 
of the enhancement of channel capacity due to entangle- 
ment and the study of multipartite entanglement. For a 
more detailed motivation the interested reader may con- 
sult [J. 

The paper is organized as follows : In Sec. |TI]we sum- 
marize the main results of [l] and then introduce the 
concept of Gauss sums and some useful related proper- 
ties. Next we provide an explicit parametrization of a 
family of equientangled bases and we prove that it inter- 
polates continuously between a product basis and a max- 
imally entangled basis, for all dimensions. We illustrate 
the behaviour of our solution with explicit examples. In 
Sec. IIIII we construct another such family using a com- 
pletely different method based on graph states, describe a 
simple extension of it to multipartite systems, and then 
illustrate its behaviour with specific examples. Finally 
in Sec. IIVI we compare the two solutions and make some 
concluding remarks. 
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A. Summary of previous work 

Let us start by summarizing the main results of [l[. 
Consider a bipartite Hilbert space T-L £g> H, where both 
Hilbert spaces have the same dimension D. The authors 
first defined an arbitrary normalized bipartite state 
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\ipo,o) = ^2 a k\k)\k). 



(1) 



k=Q 



Next for m, n = 0, 1, . . . , D — 1, they considered the col- 
lection of D 2 "shifted" states 
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|Vm,n> = X" 1 ® X 1 
D-l 

= ^a k \k(Bm)\k(Bm@n), (2) 



fc=0 



where 
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X :=J2 \k®l)(k\ 



(3) 



k=0 



is the generalized Pauli (or shift) operator and denotes 
addition modulo D. They noted that all states have the 
same value of entropy of entanglement Q given by the 
von-Neumann entropy 
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a-kl 



(4) 



k=Q 



where the logarithm is taken in base D for normaliza- 
tion reasons so that all maximally entangled states have 
entanglement equal to one regardless of D. 

Demanding the states in $2$ be orthonormal yields 
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{ak)*a k <s,m = S m> o, Vm = 0, . . . , D - 1, (5) 



fc=0 
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and the authors proved (see their Eqn. (36)) that ([5]) is 
satisfied if and only if the coefficients a./, have the form 
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(ik 



D ^ 



(6) 
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where 0j are arbitrary real parameters and lu = e 27n / D is 
the D-th root of unity. 

Therefore the authors found a family of D 2 orthonor- 
mal states, all having the same Schmidt coefficients and 
hence the same value of entanglement. To ensure it in- 
terpolates from a product basis to a maximally entangled 
basis, it is sufficient to find a set of parameters {O^^Sq 1 

for which the magnitude of ak is \a,k\ — 1/V~D for all fc. 
Then the problem is solved by defining 



D-l 



(7) 



where t € [0, 1] is a real parameter. When t = we have 
dk = 5k,o so the basis states are product states and when 
t = 1, the basis is maximally entangled by assumption. 
We also observe there is a continuous variation in between 
these two extremes as a function of t. 

Karimipour and Memarzadeh considered the existence 
of such a set {Oj}®^ 1 in arbitrary dimensions (see the 
last paragraph of Sec. V in They found particular 
solutions for D ^ 5, but did not find a general solution 
for arbitrary D. 



B. Quadratic Gauss Sums 

We now define the basic mathematical tools we will 
make use of in the rest of this section. The most impor- 
tant concept is that of a quadratic Gauss sum, defined 
below. 

Quadratic Gauss Sums. Let p, m be positive integers. 
The quadratic Gauss sum is defined as 



Generalized Quadratic Gauss Sums. Let p, m, n be 

positive integers. The generalized quadratic Gauss sum 
is defined as 
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27ri(j m-\-jn) jp 



(10) 



Finally the following reciprocity formula for general- 
ized Gauss sums holds. 

Reciprocity Formula for Generalized Quadratic 
Gauss Sums. Let p, m, n be positive integers such that 
nip =/= and mp + n is even. Then 



1 p 

jp 



^7ri(mp — n 2 )f 4mp 
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-iri(j 2 p+jn)/m 



(11) 



3=0 



The definitions of the Gauss sums (jHJ) and (|T0|) as well 
as the Landsberg-Schaar's identity J3) can be found in 
standard number theory books [3|-|5(- The reciprocity 
formula for the generalized quadratic Gauss sum is not 
as well-known, and can be found in Q. 



C. Explicit Solution 

We now show that a family of equientangled bases that 
interpolates continuously between the product basis and 
the maximally entangled basis exists for all dimensions 
D, as summarized by the following Theorem. 

Theorem 1. The collection of D 2 normalized states 

D-l 

\ip m ,n{t)) = a k (t)\k®m)\k®m®n), (12) 

m, n = 0, . . . , D — 1, 

indexed by a real parameter t G [0, 1] with 
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4=0 



(8) 
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The quadratic Gauss sums satisfy a reciprocity relation 
known as 

Landsberg-Schaar Identity. Let p, m be positive inte- 
gers. Then 
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e 27rij m/p _ 
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3=0 



-nij p/2m 



(9) 



The quadratic Gauss sums can be generalized as fol- 
lows. 



<*(*) = s E ^ 



LO = C 2 ^' D , 



3=0 



with the particular choice of 



-jo J ft] /D if D is even 
i = 1 2TTj 2 /D if D is odd, 



(13) 



(14) 



defines a family of equientangled bases that continuously 
interpolates between a product basis at t = and a max- 
imally entangled basis at t = 1 . 

That (|12p defines a family of equientangled bases that 
consists of a product basis at t = follows directly from 
the remarks of Sec. HI a^(0) = Sk,o- Next note that a 
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continuous variation of t in the interval [0, 1] corresponds 
to a continuous variation of the Schmidt coefficients of 
the states in the basis. The latter implies that no matter 
which measure one uses to quantify the entanglement, 
the measure will vary continuously with t, since any pure 
state entanglement measure depends only on the Schmidt 
coefficients of the state Q. 

The only thing left to show is that the basis states 
in Theorem [JJ are maximally entangled when t = 1, or, 
equivalently, that |<Zfe (1) | = 1/VD for all k. We prove 
this by explicitly evaluating the value of afe(l) in the 
following Lemma. 

Lemma 1. Let au(t) and {O^fSy 1 be as defined by The- 
orem [TJ Then for all k 



D=5 



afe(l) 



a 7Ti/4 



-k 2 /2 



D w 



-k- 



if D is even 
74 H^)< jf-D is odd ■ 

(15) 



Lemma [JJ implies at once that |afc(l)| = 1/yD, and 
therefore proves Theorem [TJ 

Proof, (of Lemma [JJ Note first that the expression for 
afc(l) in ([TBI with 6® defined in (|TTj) resembles the gen- 
eralized quadratic Gauss sum (fTU)) . We will use the reci- 
procity formula (jllj) to prove Lemma [JJ There are two 
cases to be considered: Even D and odd D. 

Even D. Note that one can rewrite afc(l) in (fT3|) with 
0° defined in (TH} as 



D-l 



1 Z5 - 1 

— e m(j 2 +2fej)/-D^ 



J=0 



j=0 



(16) 

Applying the reciprocity formula (fTTj) to last term in (jTSJ) 
with m = 1, n — 2k, p — D (noting that mp + n = D + 2k 
is even) yields 

fflfc (l) = ^ c -p-4fc 2 )/4Z5 = i^^^Cc-.rifc 2 /!? 



3 iri/4 



,-fe 2 /2 



since (— 1) = 1 for even D. 



(17) 



Odd 13. The proof is essentially the same as in the 
even D case, but we explicitly write it below for the sake 
of completeness. Using a similar argument we rewrite 
Ofc(l) in {13} with 9° defined in ^ as 



D-l 1 D-l 



3=0 



3=0 



(18) 



Applying again the reciprocity formula pT|) to last term 
in (|18p with rri = 2, n = 2fc,p = Z) (noting that mp + n = 




FIG. 1. The variation of |afc(i)| [dimensionless] with t for 
D = 5. Note how at t = all coefficients but one are zero, 
and how at t = 1 all coefficients are equal in magnitude to 
1 / y/5, with a continuous variation in between. The horizontal 
dashed line represents the l/v5 constant function. 



2D + 2k is even) yields 

afc (l) = _Lf (1 + e -iri(D+2fc)/a) 



V2 



7ri/4 /I ;2fc+D 

e ,-fc 2 /4 ' 1 1 

rUJ 



V2 



(19) 



where we used (— \) 2k+D = _(j2fe+£ ) 'j smce 2k + D is odd. 
This concludes the proof of Lemma [1] and implicitly of 
Theorem CD □ 



D. Examples 

In this section we present some examples that illustrate 
the behaviour of the solution we provided in Theorem [TJ 
for various dimensions. First we consider D = 5 and we 
plot the absolute values of the (t) coefficients as a func- 
tion of t in Fig. [TJ It is easy to see that indeed the basis 
interpolates between a product basis and a maximally 
entangled one in a continuous manner. We observe that 
all coefficients are non-zero for t > and we believe that 
this is probably also the case for all odd D's. 

In Fig. [3] we perform the same analysis as above, but 
now for D = 8. We observed that some coefficients vanish 
for some values of t, which seems to be true in general 
for even D. 

In Fig. [3] we plot the entropy of entanglement of the 
states in the basis as a function of t for dimensions D — 
2,3,5,8 and 100. We see how the entanglement varies 
continuously but not monotonically between and 1. 

Finally in Fig. [4] we display a parametric plot of the 
variation of the second Schmidt coefficient ai(t) in the 
complex plane as t is varied from to 1 for D = 51, so 
that the reader can get an idea of how the coefficients 
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D=8 

|a*(t)l 




FIG. 2. The variation of |afc(i)| [dimensionless] with t for D = 
8. Again note how at t = all coefficients but one are zero, 
and how at t = 1 all coefficients are equal in magnitude to 
l/v8, with a continuous variation in between. The horizontal 
dashed line represents the 1 / y/S constant function. 

Normalized 
Entropy of 
Entanglement 




FIG. 3. The normalized entropy of entanglement [dimension- 
less, the logarithm is taken base D] as a function of t for 
various dimensions. Note that the variation is not monotonic 
(except for D — 2), although for large D the oscillations tend 
to be smoothed out. 



denned in (fT3| look in general. The other coefficients ak 
look similar. 



III. CONSTRUCTION BASED ON GRAPH 
STATES 

A. Explicit solution 

We provide below another solution to the problem that 
uses qudit graph states. Again having in mind a bipartite 
Hilbert space T-L'&'H, both local spaces having dimension 



D=51 
Im[a,(t)] 




FIG. 4. Parametric plot of ai(t) [dimensionless] in the com- 
plex plane as t is varied from to 1. Note that ai(0) =0 and 
a\(l) = -J==e' rl '3 _ 2~5T', the value provided by Lemma U 
The starting point t = and the ending point t — 1 are 
marked by solid disks. 

D, we define a one-qudit state 

y/D k=0 

It is easy to see that the collection of D states 

\m):=Z m \+), m = 0,...,D-l (21) 

defines an orthonormal basis of T~L (also known as the 
Fourier basis), (m\n) — S mn , where 

Z:=5> fe |fc)(fc|, (22) 

fc=0 

with ijj = e 2m / D being the D-th root of unity. It then 
follows at once that the collection of D 2 states 

\m)\n) - (Z m ®Z n )\+)\+), m,n = 0,...,D-l (23) 

defines an orthonormal product basis of the bipartite 
Hilbert space Ht^H. 

Next we define the generalized controlled-Phase gate 

as 

D-l D-l 

CP := £ \k)(k\ ®Z k =Y. ^Ii>01 ® \k){k\ (24) 

k=0 j,k=0 

and note that CP is a unitary operator that commutes 
with Z m <g> Z n , for all m, n = 0, . . . , D - 1. The state 

|G):=CP|+)|+) = i f>^|j)|fc) (25) 

j,h=0 
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is an example of a two-qudit graph state and it is not 
hard to see that \G) is maximally entangled. Then the 
collection of D 2 states 



(Z m <gi Z n )\G) = (Z m <E> Z n )CP\+)\+) 
= CP(Z m ®Z n )\+)\+), m,n = 0,. 



,D-1 (26) 



defines an orthonormal basis of the bipartite Hilbert 
space T-L®T-L, which we call a graph basis. Since Z m ® Z n 
are local unitaries and \G) is a maximally entangled state, 
then all the other graph basis states must also be maxi- 
mally entangled. For more details about graph states of 
arbitrary dimension see [ij-tioj. 

We now have all the tools to construct a continuous in- 
terpolating family of equientangled bases, as summarized 
by the Theorem below. 

Theorem 2. The collection of D 2 normalized states 

\G m , n (t)) = (Z m ® Z")CP(t)|+)|+>, (27) 
m, n = 0, . . . , D — 1, 

indexed by a real parameter t € [0,1] where 

CP(t)=J2" 3kt \j)(3\®\k)(k\ (28) 

defines a family of equientangled bases that continuously 
interpolates between a product basis at t = and a max- 
imally entangled basis at t — 1. 

Proof. We make the crucial observation that CP(t) com- 
mutes with Z rn <£) Z n for all m, n = 0, . . . , D — 1 and all 
t e [0, 1] which implies that {\G m ^ n (t))}^l =0 defines an 
orthonormal basis since it differs from the orthonormal 
basis in (|2"3")1 only by the unitary operator CP(i). All 
states in the basis are equally entangled, and moreover, 
share the same set of Schmidt coefficients since any two 
basis states are equivalent up to local unitaries of the 
form Z m <g) Z n . 

Finally note that CP{t = 0) = I ® I and 
CP(t = 1) = CP (defined in (|23))). and therefore at 
t = the basis is product, see (j23|) . and at t = 1 the 
basis is maximally entangled, see (|2"6"]l . The operator 
CP(t) can be viewed as a controlled-Phase gate whose 
"entangling strength" can be tuned continuously. The 
Schmidt coefficients of the states in the basis vary 
continuously with t and hence the entanglement also 
varies continuously with t, regardless of which entan- 
glement measure one uses (see the remarks following 
Theorem [[]). □ 

Our construction above can be expressed in the frame- 
work described in the last two paragraphs of Sec. II of 
[H, by setting U m = Z m and V„ = Z n . 

Next we prove that the Schmidt coefficients of the basis 
states in Theorem [2] are all non-zero for any t > 0, so all 
bases consist of full Schmidt rank states whenever t > 0. 

Lemma 2. The equientangled family of bases 
{\G mi n(t))}m~n=o defined in Theorem [H consists of 
full Schmidt rank states, for any < t ^ I . 



Proof. We will show that the product of the Schmidt coef- 
ficients is always non-zero, which implies that no Schmidt 
coefficient can be zero, whenever < t ^ 1. 
Let 



D-1 



(29) 



j,k=0 



be an arbitrary normalized pure state in a bipartite 
Hilbert space H <£> H and let {A^} denote the set of 
Schmidt coefficients of satisfying J2k^k = 1; note 
that they are equal to the squares of the singular val- 
ues of the coefficient matrix f2 in (|29p . The product of 
the squares of the singular values is just the product of 
the eigenvalues of , the latter product being equal to 
det(ilfit) = | det (il) | 2 , so we conclude that 



D-1 



11 X k = |det(fi)| 2 



(30) 



k=Q 



The states \G„ l:n (t)) of Theorem [2] share the same 
set of Schmidt coefficients (they are all related by lo- 
cal unitaries) so it suffices to show that the product of 
the Schmidt coefficients is non-zero only for the state 
\G ,o(t)). Recall that 



|G 0)0 (*)) = CP(t)|+)|+) 



E 



jjkt 

~D 



-\j)\k). 



(31) 



Expressing the coefficients u)i /D as a matrix fl(t), one 
can easily see that D ■ ft(t) is a D x D Vandermonde 
matrix whose determinant is 



det [D-O(i)] = Y[(u jt 

j>k 



,kt.\ 



j>k 



kt 



,U~k)t 



1 



(see p. 29 of 11( for more details on Vandermonde matri- 
ces). For a given t, the product above is zero if and only 
if at least one term is zero, i.e. there must exist integers 
j, k, with 0^k<j^D — I, such that 



(j - k)t = nD t = n 



D 

j-k' 



(33) 



for some positive integer n ^ 0. Note that < j — k ^ 
D— I , so -D/ (j — k) > 1 and the above equation can never 
be satisfied for < t ^ I . We have therefore proved that 
det[fl(t)] £ for < t < 1, which, in the light of (|3TJ|). 
is equivalent to saying that the product of the Schmidt 
coefficients is non-zero for < t ^ 1, and this concludes 
the proof of the Lemma. □ 

For this family of equientangled bases, we do not have 
an analytic expression for the Schmidt coefficients nor 
the entropy of entanglement for general D though they 
can be easily found by numerically diagonalizing the co- 
efficient matrix £7(t)Sl(t)^. Having said that, we derived a 
simple analytic expression for the product of all Schmidt 
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c g M) = d( 11 a 



coefBcients, see (|30j) and (|32j) . which is simply related 
to an entanglement monotone called G-concurrence (first 
introduced in [l2j]) which is defined for a pure bipartite 
state (|29p in terms of its Schmidt coefficients {A*,} as 

k) =D\det(n)\V D (34) 

\fe=0 / 

where J^ fc Xk = 1. The G-concurrence is zero whenever 
at least one Schmidt coefficient is zero and is equal to one 
if and only if the state is maximally entangled. Unlike 
the entropy of entanglement, we are able to show two 
analytical results that are true for all D expressed in 
Lemmas [3] and 2] 

Lemma 3. The G-concurrence of the equientangled basis 
states, {|G mj „(i))}^^ l 1 =0 defined in Theorem\^is 



C G (t) = 



D 



D-l 



Yl [sin 2 (wrt/D)] 



(D-r)/D 



for all m, n, D. 

(35) 



Proof. Every basis state has the same G-concurrence 
since they all share the same set of Schmidt coefficients 
(recall that they differ only by local unitaries). Invoking 
the definition of G-concurrence, we have 



C G (t) = D\det(fl(t))\ 2/D = D 



D D 



det(D -n{t)) 



2/D 



D 



n 

3>k 



,U-k)t 



2/D 



D 



D-l 

n 

r=l 



D-l 



in 

r—1 

i D - 1 



} (D-r)/D ^ _ C0S ( 2 TTrt/D)} — 
~~ r=l 

1 J-J 1 2 2(D-r)/D [ sin 2 (7rrt/jD) ] ^ 

r=l 

D-l 

2/__ J ,mlP-'')/D 



,D-1 



D 



Yl [sin 2 (7rrf/L>) 



(36) 

2(£>-r) 

75 

(37) 
(38) 

(39) 

(40) 

(41) 



r=l 



where in (|36|) we used the fact that det (cM) — c D det(M) 
for a D x D arbitrary matrix M and an arbitrary con- 
stant c. The first equality in (|3"T)) follows at once from 
(1321) . whereas the second equality in (|37)) follows from a 
simple counting argument in which one replaces j — k by 
r, making sure that the different pairs (j, k), j > k, that 
give rise to the same r are counted; for a given r there 
are D — r such pairs. □ 

It turns out the G-concurrence has the following nice 
property: 



Lemma 4. The G-concurrence of the basis states 
{\G mt n(t))}m n=o * n PS]) is strictly increasing in the open 
interval t £ (0, 1), /or all dimensions D. 

Proof. We prove this by showing that the first derivative 
of the CG{t) with respect to t is strictly positive. First 
note that since Coif) > is a positive function in the 
interval t G (0, 1). This means showing 4fCa(t) > is 
equivalent to showing th at 4 ? logCc(i) > 0. The deriva- 
tive of the logarithm of (|4"T|) is 

d 2 

- log C G (t) = r(D-r)cot(Trrt/D), (42) 



where cot(-) denotes the cotangent function. We differ- 
entiate again to get 



dt 2 



log C G (*) 



27T 2 



1 



sin (jrrt/D) 



(43) 



Note that the right hand side of (|43l) is strictly negative 
whenever t > 0, which implies that the first derivative of 
the logarithm (|42|) is a strictly decreasing function of t 
and hence achieving its minimum value at t = 1, which 
is given by 



dt 



log G G (t) | i=1 



2vr 



D-l 



r(D - r) cot(vrr/L>) = 0, 



where the last equality follows from symmetry consid- 
erations (terms cancel one by one). We have shown 
^logC G (i) >0» gjC G (t) > and therefore we con- 
clude that the G-concurrence is strictly increasing for 
1 6(0,1). □ 



B. Extension to multipartite systems 

The construction presented in Theorem [2] can be eas- 
ily generalized to multipartite systems of arbitrary di- 
mension. The concept of maximally entangled states is 
not defined for three parties or more, so in this case, the 
family continuously interpolates between a product basis 
and a qudit graph basis. It is still true that for a fixed 
t, all basis states constructed this way have the same en- 
tanglement (as quantified by any entanglement measure) 
since they only differ by local unitaries. 

As a specific example, consider the tripartite GHZ 
state (| 000) + |lll))/\/2. This state is a stabilizer state 
and therefore is local-unitary equivalent [l3] to a 
graph state |G) = (CP 12 CP 2 3CP 13 )|+)i|+)2|+>3, where 
the subscripts on CP indicates which pair of qubits the 
CP gate is applied to. By varying the "strength" of the 
controllcd-Phase gate, one can now construct a family of 
equally entangled basis for the Hilbert space of 3 qubits 
that continuously interpolates between a product basis 
and the GHZ-like graph basis. This GHZ construction 
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D=5 D=8 
[Schmidt(t)] 1/2 [Schmidt©] 1/2 




FIG. 5. The square roots of the Schmidt coefficients [dimen- 
sionless] as functions of t for D — 5. Note how at t = all 
coefficients but one are zero, and how at t — I all coefficients 
are equal in magnitude to l/v5, with a continuous variation 
in between. The horizontal dashed line represents the 1/V5 
constant function. 



FIG. 6. The square roots of the Schmidt coefficients [dimen- 
sionless] as functions of t for D = 8. Again note how at t — 
all coefficients but one are zero, and how at t = 1 all coef- 
ficients are equal in magnitude to l/v8, with a continuous 
variation in between. The horizontal dashed line represents 
the 1 / v8 constant function. 



can be easily generalized to higher dimensions and also 
to n parties by using the complete graph given by 

n — l n 

iG GHZ (i)):=nn cp ^wi+)' 5 " ( 4s ) 

2—1 j>i 

where |+) is defined in ([20]) and CP(i) is defined in ([28]). 
Finally note that this construction works for any graph 
state of any dimension, and not just for GHZ-like graph 
states. Such bases with tunable entanglement may be of 
use in the study of multipartite entanglement. 

C. Examples 

In this subsection we perform a similar analysis as the 
one in Sec. Ill D[ so that one can easily compare the be- 
haviour of both solutions. 

We consider again a D = 5 example, for which we plot 
in Fig. [S] the square root of the Schmidt coefficients as 
functions of t. It is easy to see that indeed the basis 
interpolates between a product basis and a maximally 
entangled one in a continuous manner. As proven in 
Lemma [H all Schmidt coefficients are non-zero for t > 0. 

In Fig. [5] we plot the same quantities for D = 8. 
Observe how in both examples above the variation of the 
Schmidt coefficients is not oscillatory, as in the examples 
of Sec. HTDl 

In Fig. [7] we plot the entropy of entanglement of the 
basis states as a function of t for dimensions D = 2, 3, 5, 8 
and 100. We observe that the entropy of entanglement 
varies continuously and monotonically between and 1. 
It is not known if the entropy of entanglement is always 
strictly increasing for all D although we verified this by 
visual inspection for all D ^ 10. In Fig. [51 we plot the 



Normalized 
Entropy of 
Entanglement 




FIG. 7. The normalized entropy of entanglement [dimension- 
less, the logarithm is taken base D] as a function of t for var- 
ious dimensions. Note that the variation seems to be mono- 
tonically increasing for all D, a statement we did not prove. 

G-concurrence for the same dimensions. We see how the 
curves are strictly increasing and this is true for all D as 
proven in Lemma [4] 



IV. CONCLUSION 

We have solved the problem posed in [1] by provid- 
ing two families of equientangled bases for two identical 
qudits for arbitrary dimension D. The construction of 
the first solution is based on quadratic Gauss sums and 
follows along the lines of [l|, whereas the second family 
is constructed using a different method based on qudit 
graph states. 
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C G (t) 




FIG. 8. The G-concurrence [dimensionless] as function of t 
for various dimensions. The variation is strictly increasing in 
t for all D as shown in Lemma [4] 



of the states in the basis varies non-monotonically with 
t for D > 2. 

The second solution based on graph states consists en- 
tirely of full Schmidt rank states for any < t ^ 1 
that seem to have an entropy of entanglement that is 
strictly increasing as t increases. Unfortunately we did 
not find a simple analytic expression for the Schmidt co- 
efficients, but they can be computed numerically with- 
out much difficulty. We found a simple analytic expres- 
sion for another pure state entanglement measure, the 
G-concurrence, which we proved is strictly increasing as 
t increases. Finally we remark that one can extend this 
construction to equally entangled bases of more than two 
parties that interpolate continuously between a product 
basis and a graph basis even if the concept of maximally 
entangled states is not defined for more than two par- 
ties. This construction may be of interest in studying 
multipartite entanglement. 



The first solution based on quadratic Gauss sums has 
an explicit analytic expression for the Schmidt coeffi- 
cients that is easy to evaluate since they are just sums 
with D terms (see (fT5|) and (|14p). However some Schmidt 
coefficients can be zero and the entropy of entanglement 
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